Trapping ions with lasers 
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This work theoretically addresses the trapping of an ionized atom with a single valence electron 
by means of lasers, analyzing qualitatively and quantitatively the consequences of the net charge 
of the particle. In our model, the coupling between the ion and the electromagnetic field includes 
the charge monopole and the internal dipole, within a multipolar expansion of the interaction 
Hamiltonian. Specifically, we perform a Power-Zienau-Woolley transformation, taking into account 
the motion of the center of mass. The net charge produces a correction in the atomic dipole which 
is of order m e /M with m e the electron mass and M the total mass of the ion. With respect to 
neutral atoms, there is also an extra coupling to the laser field which can be approximated by that 
of the monopole located at the position of the center of mass. These additional effects, however, are 
shown to be very small compared to the dominant dipolar trapping term. 



I. INTRODUCTION 



Charged particles are routinely trapped in several labo- 
ratories worldwide [lj,|2||. Confinement is usually achieved 
by means of radiofrequency (rf) traps or by combining 
static electric and magnetic fields [l| ■ These technologies 
have reached such a level of precision that they are now _at 
the basis of mass spectrometers [j| , optical clocks 
and prototypes of quantum information processors 

The creation of optic al p otentials for single ions by 
means of laser beams d, [l(f has been discussed with the 
purpose of employing established techniques routinely 
used to trap neutral atoms [ll|. Recently, in Ref. [l2[ 
the trapping of a single ion in a dipole trap was reported. 
This alternative approach to confine ions optically is ex- 
pected, among other aims, to allow for the study and 
realization of controlled dynamics coupling atoms and 
ions (l3| . thereby overcoming limitations which could 
arise from the field of ion traps (T^4l7| . However, in com- 
parison to neutral atoms, ions have a net charge which 
additionally couples to the electromagnetic field. Thus, 
one should evaluate the impact of the oscillation of the 
ion in the optical potential, and the direct coupling of 
the charge to the time- varying electromagnetic field. In 
common rf traps, the latter leads to a direct drive of the 
ion's motion at the rf frequency - an effect called mi- 
cromotion [l| that can severely affect experiments that 
require not only low temperatures but also small resid- 
ual kinetic energies. Even though the physical principles 
of radiofrequency traps and optical traps are similar, the 
frequency of the electromagnetic field differs by eight or- 
ders of magnitude and therefore leads to dramatically 
different weights for the amplitudes of, in principle, sim- 
ilar effects. 

Starting from this motivation, in this work we ana- 
lyze theoretically the efficiency of laser trapping of an 



alkali-earth metal ion, namely, an ionized atom with a 
single valence electron. We start from the Hamiltonian 
in Coulomb gauge, treating the electromagnetic field in 
second quantization and the atom in the regime in which 
relativistic effects are perturbative corrections. We then 
perform the Power-Zienau-Woolley transformation, tak- 
ing into account the quantized center-of-mass motion. 
This allows us to identify a multipole expansion, in which 
the dominant term is responsible for dipolar trapping. 
We analyze the effect and estimate the order of magni- 
tude of the corrections, specially focusing on the terms 
which originate from the charge monopole. In the spirit 
of [HI , reporting a study on the stability of a single ion 
in a Paul trap, in our work we perform a stability anal- 
ysis of an ion in a dipole trap, where an effect similar to 
micromotion emerges from the coupling of the ion charge 
with the fast-oscillating laser field, and we briefly discuss 
the impact on potential applications. For definitcness, we 
shall consider the setup drawn in Fig. [TJ but our results 
remain valid for ions trapped in optical lattices. 



FIG. 1: An ion is confined in a dipolar trap formed by a 
focused red-detuned laser beam, with wave vector kL- The 
particle couples with the external fields both via the charge 
monopole Q and the optical dipole d. The first coupling gives 
rise to an effect which is similar to micromotion in an rf trap; 
the latter is responsible for the optical trapping. An electro- 
static field Eo can be used to prevent the escape of the particle 
along the propagation direction of the laser, as realized in [T3 |. 
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This paper is organized as follows: In Sec. [IT] the ba- 
sic physical concepts underlying the dipolar trapping of 
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particles are briefly reviewed. In Sec. lIIII the Hamiltonian 
in Coulomb gauge is introduced and the Power-Zienau- 
Woolley transformation is performed. The various mul- 
tipolar coupling terms emerging from the transformation 
are discussed in Sec. IIV1 while concluding remarks are 
presented in Sec. [V] 



II. DIPOLAR TRAPPING 

We first briefly summarize the concepts at the basis of 
dipolar trapping of atoms (we refer the reader to for 
a review). We assume an alkali-metal atom. For suffi- 
ciently weak laser fields, the coupling between the atom 
and the laser field only involves the valence electron, and 
can be described by the interaction Hamiltonian in the 
electric dipole approximation, which reads 



H dip = -d ■ E(R) 



(1) 



where d is the atomic dipole and E(R) is the electric field 
corresponding to the trapping laser, evaluated at the po- 
sition R of the atomic center of mass. In the following 
we will assume that the atomic motion has been brought 
to energies corresponding to temperatures of the order of 
the millikelvin by means of laser cooling [l9|, [2(| ■ In this 
regime, the characteristic time scale for the evolution of 
the atomic center of mass is usually much longer than 
the time scale in which the electronic degrees of freedom 
evolve towards the internal steady state. One can then 
evaluate the mechanical effects on the atom by means of 
an adiabatic treatment, namely, assuming that the inter- 
nal state is the stationary state corresponding to the field 
at position R, and which results from the competition of 
the dipolar excitation due to the laser and the relaxation 
due to spontaneous decay [l9j . 

We consider thus a laser which drives a transition be- 
tween a ground state \g) and an excited state |e), forming 
an optical dipole transition at frequency uj eg . The energy 
diagram is sketched in Fig. [3] The laser is a classical field, 
more specifically, a travelling wave with frequency ul, 
wave vector Ul, polarization £l, and amplitude El(R), 
corresponding to a tightly focussed beam. From the so- 
lution for the internal steady state of the atom as a func- 
tion of R, one can find the mechanical force governing the 
atomic center-of-mass motion and which reads [2lL |22j : 



<F) = -Mvin(l + S ) 



ht 



l 



(2) 



(after averaging over the optical oscillation period). Here, 
s = s(R) is the saturation parameter, which depends on 
the position through the spatial dependence of the Rabi 
frequency, and reads 



with r the natural linewidth of the dipolar transition, 
S = lul — Lo eg the detuning of the laser frequency from 
resonance, and 



fi(fl) = \ (g\d\e) ■ e L E L (R) 



(4) 



the Rabi frequency. 

The first term in Eq. @ describes a dispersive force 
which emerges from the intensity gradient, namely, the 
so-called dipolar force. The second term is the radiation- 
pressure force, which is related to dissipative processes. 
The strength of the dipolar force can exceed by orders 
of magnitude the radiation pressure when the atom is 
weakly driven, namely, when the saturation parameter 
satisfies s < 1, and when the laser is far detuned from 
the transition, such that \S\ 3> T/2 (provided that the 
field is still far-off resonance from other electronic states) . 
In this limit, the mechanical effects due to the radiation 
pressure can be neglected and the effect of the dipolar 
coupling can be described by the dipole force. This leads 
to the definition of an effective conservative potential for 
the center of mass, which is proportional to the laser 
intensity and takes the form 



V eS (R) - ^s(R) 



(5) 



(where the approximation s <C 1 has been used). Since 
s(R) is positive, the sign of the potential is solely deter- 
mined by the detuning <5. When <5 < 0, the maxima of 
the field intensity are minima of the effective potential, 
and the atom can be trapped in high-intensity regions. 
In a dressed-atom picture, the effective potential can be 
identified with the position-dependent energy shift of the 
atomic levels, as displayed in Fig. [5] 



Iff) 



Iff') 




FIG. 2: Left panel: Optical trapping is achieved by means 
of a laser that drives a dipolar transition between a ground 
state \g) and an excited state \e). The radiation is detuned 
by an amount S = uil — w e9 with respect to the atomic tran- 
sition. Right panel: In the dressed atom picture, correspond- 
ing to the eigenstates of the atomic Hamiltonian including 
the coupling to the classical laser field, the levels have energy 
shifts that depend on the Rabi frequency Q. This gives rise 
to an effective potential associated with the spatial variation 
of the laser intensity, which is proportional to Q 2 and is here 
reported along a direction x transverse to the propagation 
direction of the laser. 



(3) 



Radiation pressure can slightly shift the equilibrium 
position with respect to the minimum of the effective 
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potential. Moreover, scattering of photons by the atom 
leads to heating of the atomic motion by an amount of 



the Coulomb (transverse) gauge takes the form 



the order of the recoil energy E T , 



(hk L ) 2 /2M per 



scattering event. The rate at which photons are scattered 



is given by [11 



s h ' 



(6) 



where the formula is valid at low saturation. For |<J| 3> 
r/2, a time scale of the dynamics can then be identified, 
in which atomic motion is essentially determined by the 
conservative potential V c g , while incoherent scattering of 
laser photons can be neglected. Further orders of mag- 
nitude in the ratio between potential depth and heating 
may be gained by using a blue-detuned standing wave 
laser, such that the atoms are confined in the nodes of 
the field, as in dispersive optical lattices 23 1. 

In the following we analyze how these dynamics are 
modified when the atom has a net charge which also cou- 
ples to the rapidly oscillating electromagnetic field. 



III. ION-PHOTON INTERACTIONS 

The system we shall consider is an ionized alkali-earth 
metal atom, with mass M and total charge Q. We treat 
the ion as a hydrogen-like system, since the radiation 
which illuminates it is optical and at sufficiently low in- 
tensity, such that the induced electronic transitions in- 
volve only the valence electron. Center-of-mass and elec- 
tronic excitations will be treated systematically in the 
non-relativistic limit. We will focus on the mechanical 
effects of light on the center-of-mass motion, considering 
both the coupling of the field with the charge monopole 
and with the atomic dipole. 

For the sake of simplicity we shall take a spinless nu- 
cleus as in the experiment reported in [l2j | , therefore not 
considering any hyperfine structure 1 . 



(i). Hamiltonian in Coulomb gauge 

We assume a hydrogen-like ion with a core of mass 
m n and charge q n , and we denote by m e and q e the mass 
and charge of the valence electron respectively. The total 
charge is thus Q = q n + q e , and the total mass is M — 
m n +m e . In the non-relativistic regime the Hamiltonian 
for the atom and the quantum electromagnetic field in 



1 In general, due to the large mass difference between electron and 
nucleus, the coupling of the magnetic field to the nuclear spin is 
much weaker than that between the magnetic field and the spin 
of the electron, so that it gives rise to corrections of higher order 
than the ones reported here. 



j=n,e 

Vboul + H r ci + flj-ad (7) 



where fj and pj (j = n, e) are conjugate position and mo- 
mentum, A is the vector potential in the Coulomb gauge, 
Vboul contains the Coulomb interaction between core and 
electron, $ represents an external electrostatic potential, 
.ffrei gives the relativistic corrections to the atomic Hamil- 
tonian, and -ffrad is the energy of the electromagnetic ra- 
diation in vacuum, 



ad 



(8) 



Here the subscript m runs over the modes of the elec- 
tromagnetic field, characterized by the frequency uj m , 
wave vector k m , and polarization given by the unit vec- 
tor e m _L k m . The operators a m and annihilate and 
create, respectively, a photon in mode m and satisfy the 



commutation relations [<z„ 



In Eq. 



have taken a quantization volume in a cubic box of side 
L with periodic boundary conditions, and we have set 
to zero the energy of the vacuum state. In terms of the 
creation and annihilation operators, the vector potential 
is given by: 



A{x) 



E 



(2e uj m L 3 ) 



ik m -x 



H.c. 



(9) 



In order to perform a multipolar expansion, it will be 
convenient to switch to a description in terms of center- 
of-mass and relative coordinates, that we shall denote by 
R, r respectively; the conjugate momenta will be denoted 
by P and p, and we will use jj, for the reduced mass. 



(ii). Power-Zienau-Woolley transformation 

In order to study the dynamics induced on the ion by 
weak fields in the optical regime, we shall use the Dirac- 
Heisenberg line gauge [24], . We shall then perform a 
multipolar expansion in powers of the small parameter 
kr, where r is the size of the atomic bound state and 
k = 2tt/X with A the wavelength of the modes under 
consideration, of the order of the laser wavelength (for 
the experiment in kr ~ 10~ 3 ). 

Starting from Eq. ([7]), one obtains the Hamiltonian in 
the new representation by means of the Power-Zienau- 
Woolley unitary transformation [24], [26[. Although the 
two representations are equivalent, the minimal coupling 
Hamiltonian ([7} is less suitable to express the problem of 
the atom interacting with an external field in a truncated 
basis of eigenstates of the atomic Hamiltonian [2^, [28[ . 
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In particular, in the Dirac-Heisenberg line gauge the mo- 
mentum p corresponds, to zero order in the expansion in 
kr, to the mass times the velocity for the relative motion. 

We note that this change of representation is com- 
monly used in problems which restrict to the dipolar ap- 
proximation, and/or the case of neutral particles, and/or 
assume that the center of mass is fixed [23 - 1261 . [29| . In 
this paper we deal with a charged particle taking into 
account the motion of the center of mass and expanding 
in powers of kr beyond the dipolar approximation, there- 
fore considering a more general framework than the one 
usually encountered in the literature. 

The Power-Zienau-Woolley transformation is defined 
by the unitary operator T, which reads: 



where 



S 



-iS/h 



dVP (x) ■ A K (x). 



(10) 



(11) 



and the integration is over all space. The polarization 
operator Pn depends on the atomic variables, and takes 
the form 1261: 



j=n,e 



R) 



dX S 



(12) 

while A < represents a vector potential obtained from A 
in Eq. Q by introducing a high-frequency cutoff in the 
sum over the modes. This cutoff has to be chosen suf- 
ficiently high so that the coupling between the particle 
and the modes above the cutoff is far off-resonant and 
thus irrelevant for the processes under study (apart from 
mass renormalization effects [30j, which are supposed to 
be already taken into account). 

We now briefly discuss how the operators transform 
under the unitary transformation defined in Eqs. (|10p~ 
(TT2"j) . The position operators are invariant under the uni- 
tary transformation defined by T, since they commute 
with S. The momenta, on the other hand, are trans- 
formed according to: 



Pj 



Tp 3 T^ 



Pj 



(13) 



where the gradient is taken with respect to the coordi- 
nates of particle j. This shows that the direct coupling 
between the momenta and the vector potential cannot be 
totally eliminated by this transformation. Nevertheless, 
to zero order in kr it is possible to suppress the coupling 
between the relative momentum p and the field. In con- 
trast with the case of a neutral particle, however, in the 
transformed Hamiltonian the center-of-mass momentum 
remains coupled to the vector potential A(R). 

The form in which the transformation acts on the op- 
erators associated to the different fields can be obtained 
from the commutation relations: A(x) commutes with 



A(x') and with B(x'), so that these operators are invari- 
ant under the transformation, while 

lA J (x),E k (x')] = -^5j- k (x-x') (14) 

where Sj,. is the transverse delta function 261. This im- 

JK, t J 

plies that the electric field operator is transformed as 
follows: 



E -> E' = TET^ = E- — P <. 



(15) 



Here E' is the electric field in the new representation, 
and Pq- k is the polarization of the system of charges as 
given by (|12[) . now projected onto its transverse part and 
including only Fourier modes with wavelengths below the 
cutoff. The operator E in Eq. (fl"5"|) can be identified with 
the electric displacement field (divided by eo) [Ifj]. We 
shall keep using the notation E to refer to this field, and 
E' for the electric field. 



(iii). Hamiltonian in Dirac-Heisenberg line gauge 

After applying the Power-Zienau-Woolley transforma- 
tion, the Hamiltonian W = THT^ reads 

H = H ion + iJ ra d + ^coupling + #rel- (16) 

Here -ff ra d has the same form as in Eq. flS]) and H' rcX con- 
tains the relativistic corrections in the new representa- 
tion, while 

P 2 P 2 

Hion — + -g - + Vboul + E dip (17) 

depends on the ion degrees of freedom, with E^p the 
dipolar self-energy: 

1 

which is a function of the relative coordinate r. The re- 
maining Hamiltonian term couples atomic and field de- 
grees of freedom: 



£di P = — / dV{P^f, 



(18) 



^coupling = £ Qj®(rj) + H tvs 



]=n,e 



where -fftrap reads 



#trap = ~ / dVE < ■ Pot 



P-Kqm K? cl Kl M 

M + 2/x + 2M ' [ } 



(20) 



while the other terms in Eq. (|T9| involve operators de- 
fined by: 

K G u = K e + K n , 

m n K e - m e K n 



K 1X 



M 



5 



with 



Kj = VjS-qjAii 



j = n,e. 



(21) 



We shall treat these operators in a multipolar expan- 
sion, namely, in powers of kr around the position R. The 
term H tla p defined in (|2"0|) is the one at the basis of the 
dipole force. To lowest order in kr, it corresponds to the 
dipolar Hamiltonian (TTJ), i? trap = — d ■ E + 0(kr), where 
the dipole operator takes now the form 

d= -q cS r (22) 

with q e f[ an effective charge for the dipolar coupling, 



q c s = \q e \ + 



( 

M 



(23) 



The net charge thus introduces a very small correction: 
For an atomic mass of the order of 10 proton masses, 
the difference between q e s and \q e \ is of order 10~ 4 . The 
other operators read in the multipole expansion: 



K, 



(0) 
CM 

(1) 
CM 



= -QA 



-q e fi r x 



B, 



and 



(o) _ 



0, 



K 



(2) 



m n ~m e _ s 
|(r.V^)(rxB), 



(24) 
(25) 

(26) 
(27) 
(28) 



where the superindices indicate the order in kr. In 
Eqs. ([24]) - ([28]) . all fields and field derivatives are assumed 
to be evaluated at the center-of-mass position R, and the 
notation stands for a gradient with respect to R. We 
have explicitly reported the terms up to (i) first order in 
kr and in m e /M, and (ii) second order in kr and zero 
order in m e /M, as these two small parameters are com- 
parable. 



the Hamiltonian for the electron including the coup ling 
with the laser field reads in the non-relativistic limit [3lL 



H, 



elect 



(p e - q e A) 2 q e h _ - 

h Vboul + <?e* - a ■ B - 

2m e 2m e 

(p e -q e A) 4 q e h r - 

+d ■ [E T X (p e - q e A) - (p e - q e A) X S T ]} + 
q 2 e h 2 



8m3 c 4 

eh 
8m^c 2 



{El - c 2 B 2 ) + 
{{p e -q e A) 2 ,a-B} 



(30) 



where Et denotes the total electric field, including the 
Coulomb interaction with the core, the external electro- 
static field, and the electric field of the laser. In the 
previous expression, all fields are evaluated at the elec- 
tron position, and { , } denotes an anticommutator. The 
expansion in Eq. (|3T)|) is in powers of m" 1 , which can 
be interpreted as an expansion in terms of the different 
momenta compared with m e c. The first three terms in 
(|30|) were already contained in the non-relativistic Hamil- 
tonian (J?}, while the others give the lowest relativistic 
corrections. 

The Power-Zienau-Woolley unitary transformation 
acts on the operators in the expression above as explained 



in subsection III (ii) in particular, 



Pe-qeA ->• p e +K t 



1 



, nie ^ „ nip . 
P+ M P ~ Q W A+ 2 dxB 



(31) 



where a multipolar expansion has again been performed 
for the fields, which are then evaluated at the position 
of the center of mass. The electric field must also be 
transformed according to eq. (jl5[) . In this way the extra 
couplings between the ion and the electromagnetic field 
that are contained in H' rel can be found. The effect of 



these terms shall be discussed in Subsection IV (v 



(iv). Relativistic corrections 

In the following we shall take into account the relativis- 
tic corrections to the Hamiltonian ([7]), and analyze how 
the Power-Zienau-Woolley transformation acts on them. 
In the absence of external fields, the fine structure terms 
are given by: 



Ht s — 



1 



2 m 



2„2 



Pc 



4m„ 



VL - _ h 2 

IS^l . §e + -—AV( 



Coul 



(29) 



We are interested in finding the corrections to these terms 
that occur when external fields are included. In order to 
do so, we start from the Dirac equation with a fixed core 
and apply a Foldy-Wouthuysen transformation, so that 



IV. EFFICIENCY OF DIPOLAR TRAPPING OF 
CHARGED PARTICLES 

We shall now apply the formalism introduced above in 
order to determine the efficiency of trapping an alkali- 
earth metal ion by means of lasers. 

We first remark that the dipolar interaction in Eq. (|20l) , 
in the electric dipole approximation (i.e., in zero order in 
the expansion in kr), is responsible for the optical trap- 
ping of the ion, according to the physical processes sum- 
marized in Sec.|TTJ As already mentioned, the presence of 
a total charge produces a small correction to the effective 
dipole according to formula (j2"3"]l . and hence modifies the 
potential depth and motional frequencies, as the Rabi fre- 
quency n in Eq. Q now scales with q e g. The correction 
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is of the order of the ratio m e /M between the electronic 
mass and the total mass of the ion. 

With respect to the additional electrostatic field Eq, 
the corresponding coupling with the ion is dominated by 
the charge monopole, since the typical length scale of 
this potential is several orders of magnitude larger than 
the atomic size, and the field intensities involved are too 
weak to induce significant polarization effects. 

The following subsections will be devoted to estimate 
the relative contributions of the different corrections to 
the basic dipolar trapping, and point out their similar- 
ities and differences in comparison to effects well stud- 
ied for radiofrequency traps. In order to do so, we shall 
make reference to the values of the parameters in the 



experiment reported in [12j. There, 



21 



M. 



g T ion was 



optically trapped by a laser with wavelength A = 280 
nm. The laser was a tighly focussed Gaussian beam with 
a waist radius of 7 /xm, which was circularly polarized 
and red-detuned up to \8\ — 2tt x 300 GHz with respect 
to the S1/2 O P3/2 transition, with a natural linewidth 
r ~ 2n x 40 MHz. The optical depth, given by the abso- 
lute value I V c g I of the effective potential at the minimum 
of the trap, was Uq < fog X 50 mK. The trapping fre- 
quencies at the bottom of the trap were w r adiai — 27r x 
200 kHz and w ax i a i ~ 2tt x 2 kHz for the purely optical 
setup. An electrostatic quadrupole field, with axis at 45° 
with respect to the trapping laser, provided additional 
trapping with a frequency of 2n x 45 kHz. Typical fre- 
quencies for the oscillatory motion of the atomic center 
of mass were then of the order of wo ~ 2tt x 100 kHz, 
the precise value depending on the direction considered. 
From these parameters one can identify a hierarchy of 
frequency scales, displayed in Table HJ The results for 
the orders of magnitude of the several effects analyzed in 
the following subsections are summarized in Table |TTJ 



Effect 


8E/U0 


Subsection 


Q-dependent correction to d 


1(T 4 


|III(iii)| 


Coupling in relativistic corrections 


10~ 5 


|IV(v)| 


Higher-order multipole terms 


i(r 6 


IV (iv) 


Monopole coupling —QP ■ A/M 


lO" 8 


IV (i) 



TABLE II: Orders of magnitude estimated for the most im- 
portant corrections to the dipolar trapping Hamiltonian. We 
refer to the subsections indicated for further details. Much 
smaller effects have been found as a result of the monopole 
coupling with blackbody radiation (which heats the motion 
at a rate F' w 10~ 7 Hz, as shown in Subsection IV (iii) I, and 
of the time-dependence of the optical potential (resulting in 
a micromotion with amplitude 10 -20 times smaller than the 
secular motion of the center of mass, as seen in Subsection 



IV (ii) > 



of-mass momentum P. To first order in kr, it takes the 
form: 



Ti' RcM 



(32) 



The first term is in zero order in kr and describes the cou- 
pling between a charge monopole and the electromagnetic 
field, while the second term is in first order and involves 
both the center-of-mass and the dipolar degrees of free- 
dom. 

In order to estimate the magnitude of the zero-order 
term we take the following approach: We focus on the 
center-of-mass motion and make a harmonic approxima- 
tion for the effective potential V e g(R), namely, we make 
a second-order Taylor expansion about the equilibrium 
position of the center of mass. This system has then an 
effective Hamiltonian: 



Frequency 


Magnitude (2n x Hz) 


wo, E TCC /li 


(0.5 - 2) x 10 5 


r 


0.4 x 10* 


Uo/h 


10 9 




up to 0.3 x 10 11 




0.3 x 10 12 


ML, 0J e g 


10 15 



(P - QA) 2 Muj 2 R 2 



TABLE I: Orders of magnitude for the different frequ ency 
scales involved in the problem, taking parameters from [12 ] . 
The symbols have been defined in the text. 



(i). Coupling of the electric field to the total charge 

Our primary interest is the study of effects due to the 
net charge of the trapped particle. Thus, we shall con- 
sider first the term in -^coupling m which the monopole 
appears to the lowest order in the multipolar expansion, 
namely, the term that couples the field with the center- 



2M 



(33) 



where A is the vector potential for the laser field, and 
where the motional frequency ujq is given, for the purely 
optical trapping, by the laser intensity and the waist of 
the beam at the focus. For the sake of simplicity we have 
assumed the same trapping frequency for all spatial direc- 
tions, since we only wish to make an order-of-magnitude 
estimation. To this effect, we shall take ujq of about ten 
orders of magnitude smaller than the frequency of the 
laser field. In the previous effective Hamiltonian (|3"3")l 
we have also included the term {QA) 2 /(2M), the zero- 
order contribution to the terms in -f/coupiing which are 
quadratic in the field operators. Equation (1331) describes 
a harmonic oscillator driven by a field with frequency 
ljl 3> ujQ. The effect of the coupling is then to induce a 
micromotion-like oscillation with very small amplitude at 
the frequency of the laser. Typical energies associated to 
this oscillation are of the order of (QA) 2 /(2M), which is 
about eight orders of magnitude smaller than the optical 
depth, and four orders smaller than the motional energy 
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scale given by Hujq (we note that in conventional rf-traps 
the energies associated to secular motion and micromo- 
tion are normally of the same order of magnitude). 

The following term in the coupling to the total mo- 
mentum (l32l) is of order P/ (Mc) < 10~ 8 when compared 
to the dipolar trapping term (JTJ). Besides, since B and 
E are out of phase, this term is proportional to the com- 
ponent of the dipole oscillation that is out of phase with 
respect to the electric field, and which is smaller by a 
factor r/<5 ~ 10~ 4 with respect to the in-phase compo- 
nent of the dipole oscillation, associated with the optical 
trapping. Taking both aspects into account, this term 
can then be neglected. 



(ii). Time-dependence of the effective potential 



The previous considerations have taken only the op- 
tical potential into account. In a setup as the one in 
Fig. [U in which a static quadrupole field E$ is also in- 
cluded, this field will modify the value of the parameter 
a in the corresponding Mathieu equation (because the 
static trapping potential along the propagation direction 
of the laser is obtained at the cost of a repulsive trans- 
verse force). However, as long as this electrostatic force 
is weak enough compared to the dipole force, as was the 
case in the experiment reported in [12| , the conclusions in 
the previous paragraph remain valid. Finally, it is worth 
noting that the micromotion frequency resulting from the 
time dependence of the potential ([Ml) is of about 2cj£, 
and thus far-off resonance from the force exerted by the 
laser on the net charge (which was analyzed in subsection 
IV (i) and has frequency lul)- 



The effective potential V c g(R) for the center of mass 
results from averaging the optical potential over the fast 
oscillation period of the laser field. Indeed, the time- 
dependent optical potential experienced by the center of 
mass can be written as: 

^optical [R,t) = 2V cS {R) cos 2 \uj L t) = 

= V eS (R)[l + cos(2u L t)}, (34) 

as the dipolar force is proportional to the product of the 
atomic dipole and the electric field, each of them oscil- 
lating at the optical frequency lo^. One might wonder 
whether by including the time dependence of this po- 
tential some other relevant micromotion-like effect could 
appear. 

Once again we can estimate the importance of this ef- 
fect by considering a harmonic approximation with trap- 
ping frequency ujq for V e g (R). Under this approximation, 
the time-dependent optical potential (|34|) is of the same 
form as the one for the motion of an ion in each of the 
transverse directions in a Paul trap: 



Vp au l(x) 



Mx 2 U) 



RF I 



2qcos{uj RF t)} 1 (35) 



which leads to a Mathieu equation for the coordinate x. 
By comparing the two previous potentials one finds that 
the parameters a and q for the dipole trap are both of 
the order of (uiq/ujl) 2 ~ 10~ 20 . This quantity determines 
the order of magnitude for the ratio of the high-frequency 
components with respect to the secular motion in the so- 
lution of the Mathieu equation (and the parameters are 
in the stability region since, for the purely optical poten- 
tial, a > 0) [33j. The micromotion amplitude associated 
with the time dependence of the optical potential is thus 
negligible. We note that because in this case both param- 
eters a and q scale as (wo/wl) 2 , the kinetic contribution 
coming from the micromotion can also be safely ignored 
(as opposed to the results discussed in [§3] for a trapped 
ion with a = and q — > as u>rf oo, since in the case 
a = the micromotion's kinetic energy has a different 
scaling) . 



(iii). Coupling between the center-of-mass motion 
and the blackbody radiation 

The coupling between the net charge and the external 
field could become important if field modes with frequen- 
cies similar to the motional frequencies of the center of 
mass are taken into account. To study this effect, which 
could give rise to additional heating mechanisms, we shall 
consider once more the motion of the center of mass in 
harmonic approximation for the effective potential, and 
analyze the coupling between this charged harmonic os- 
cillator and the continuum of modes of the electromag- 
netic field. This problem can be modelled by the Hamil- 
tonian H c g + i? ra d with H c g introduced in Eq. (f33|) . and 
H la ,d defined in ([5]). Following [l2|], we shall be inter- 
ested in motional trapping frequencies ranging from 10 
kHz and 1MHz. For these values, the wavelength of the 
field modes which are close to resonance is many orders 
of magnitude larger than the typical size of the atomic 
motion, and then the field can be treated as spatially 
constant. 

We wish to estimate the rates at which energy is ex- 
changed between the harmonic oscillator and the contin- 
uum of modes, assuming that the electromagnetic field 
is in a thermal state. This calculation is carried out (for 
instance) in [22| , where the rate of heating of the motion 
due to the interaction with the thermal bath is shown 
to be given by a quantity T' associated to stimulated 
transitions, and thus dependent on (ti(wo)}, the mean 
population of the modes of the continuum which are res- 
onant with the motion (the rate of energy exchange also 
depends on the energy of the charged particle, but this 
is negligigle in our problem) . For ujq of the order of 2ir x 
100 kHz, the quantity HuJo/kB corresponds to a tempera- 
ture of a few fiK. The mean population of the modes res- 
onant with ujq for a radiation bath at room temperature 
T is thus (n{u>o)) w (k B T)/(huj ) w 10 s . The resulting 
heating rate is V « 10~ 7 Hz, and the associated time 
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scale is of a few months. 



which is of order of (kr) 2 compared to the zero-order 
term — d ■ E. 



(iv). Higher-order corrections in the coupling 
between the laser and the relative motion 



2. Coupling of the field with the relative momentum 



As seen in the previous subsections, the monopole cou- 
pling between the center of mass and the electromagnetic 
field gives rise to very small effects. We shall then study 
the dominant higher-order corrections to the dipolar in- 
teraction which couples the electronic degrees of freedom 
and the field, and conclude that the most important cor- 
rections are also present when considering the dipolar 
coupling of a neutral atom, so that no relevant additional 
effects that depend on the total charge of the particle ap- 
pear. 



1. Multipole expansion of Htr&p 

We first focus on Hamiltonian ([20]) . which at lowest or- 
der in the expansion in kr (electric dipole approximation) 
determines the optical potential, and discuss now the 
higher-order corrections in powers of kr, taking the low- 
est order in the expansion in powers of the ratio m e /M. 
The electric quadrupole correction reads 



H, 



(i) 

trap 



<7e 

2 



Vfl)(r-£) 



(36) 



This term is far-off resonance: The operator r 2 does 
not connect the levels involved in the dipolar transi- 
tion driven by E, because they have different parity; the 
non-vanishing matrix elements in H^ p will thus oscil- 
late too fast to contribute significantly to the evolution. 
Indeed, the effect of the coupling ([36j) can be estimated 
by means of time-dependent perturbation theory. First- 
order processes will give rise to transitions in the dressed- 
state basis with probability amplitudes of the order of 
(kr)n/ojr, ~ 1CP 8 , and second order processes will have 
associated rates that scale like (kr) 2 fl 2 /ojt,. These con- 
tributions are then smaller than the following order in the 
expansion of i/trap in powers of kr, which we denote by 
fltrap- The dominant correction to the dipolar coupling 
is then given by the electric octupole of (|2"0"|) . 



H, 



(2) 
trap 



(37) 



2 One can also check whether, even if the heating rate is very low, 
the energy shift of the levels of the harmonic oscillator due to the 
coupling with the electromagnetic field can be significant. Along 
the lines of [30ll , in order to avoid divergences a renormalization 
term can be included accounting for the energy shift of a free 
particle, and then a relativistic cutoff of order Mc 2 is introduced 
in the coupling. The result for the energy shift is of about twenty 
orders of magnitude smaller than the zero-point energy fiwo of 
the harmonic approximation for the motion. 



We next analyze the term in H coup \i ng , Eq. (|19p . which 
describes the coupling of the external field with the 
relative momentum, p ■ K Te \/ fj,, with K Te \ reported in 
Eqs. (|2"6")) - (|2"51) . In lowest order in kr, and neglecting 
terms of order (kr)(m e /M) 2 , it reads 



qe 
' 2 



1 



1 

m r , 



l-B. 



(38) 



where I = f x p is the orbital angular momentum of the 
electron. This term describes the coupling between the 
magnetic dipole, induced by the orbital motion of the 
valence electron, and the magnetic field, and it connects 
bound states with equal parity. It is thus far-off reso- 
nance from the transitions it couples, and hence negligi- 
ble. In presence of static magnetic fields, it may give rise 
to energy shifts whose magnitude depends on the angular 
momentum. 

The following terms in powers of kr have resonant con- 
tributions which must be estimated. The dominant part 
(i.e., to leading order in m e /M) is given by the expres- 
sion: 



?(2) 



K, 



P- 



2fi 



3m, 



■(r.f R )(l.B) 



(39) 



Similarly to Q37p. this term is of order (fcr) 2 when com- 
pared to the dipolar coupling, and to this order it does 
not depend on the total charge Q. 



3. Effect of the previous corrections 

The terms so far discussed constitute corrections to the 
coupling between electric field and atomic dipole, so that 
the Rabi frequency defined in Q is modified according 
to: 



Sl(R) = q cS (g\r-E\e) 



Qes(g\r- E\e)- 



(r-V R )(r-E) + —(l-B) 



\e) 



, (40) 



which is reported in rotating- wave approximation and up 
to second order in kr and first order in (kr)(m e /M). As 
a consequence, the resulting dipolar potential, that de- 
pends on the Rabi frequency as shown in formulas ([5])-([5|) 



can be written as V c s = + V^' , where the super 
script indicates the order in kr, with 



r(2) 



Kff } « \q eS (g\r-E\e)\ 



(41) 
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the dipolar potential in the electric dipole approximation, 
and 



(v). Spin coupling and corrections to the fine 
structure terms 



V$> cx fit, 



{<■ 



g\f-E\e) 



(9\[(r 
2 



V n ) 2 (r -E)+ 



+ —(f-^n)(f'S)]\e)} (42) 



the first correction in the multipole expansion. Its or- 
der of magnitude can be estimated by considering that 
it can be rewritten in terms of second derivatives of the 
electric field, and hence in terms of the typical length 
scale over which the electric field changes. In particular, 
if the electric field is sufficiently smooth in the trans- 
verse plane, then the dominant corrections are due to 
the derivatives along the axial direction, which are pro- 

(2) 

portional to k^. The resulting terms in V~J will then be 

proportional to the light intensity, just as VLgP, and will 
only provide a global factor giving a relative difference of 
order (kr) 2 ~ I0 -6 in the values of the potential depth 
Uq and the motional frequencies with respect to the ones 
obtained from the lowest-order term in the potential. If 
instead the transverse derivatives are not negligible, then 
the corrections arising from them may slightly modify the 
shape of the effective potential (an effect which can be at 
most of the same order of magnitude as the one related 
to the axial derivatives). 



4- Terms quadratic in the field operators 

Finally, there are some additional terms in i? C oupiing 
which are quadratic in the field. The zero-order in kr, 
of the form (QA) 2 /(2M), couples only to the center of 

The 



mass and has been discussed in Subsection IV (i) 



next two terms in the expansion in powers of kr are (to 
lowest order in m e /M) given by: 



Q-(dx B) ■ A+ —{dx B) 2 . (43) 



To conclude our analysis, we consider the couplings 
coming from the relativistic terms in H' rcl , obtained in 



Subsection III (iv) Typical experiments include weak 
magnetic fields, used to set a preferred quantization axis 
for the angular momentum. The magnetic field of the 
laser could in principle induce transitions between dif- 
ferent states and thus interfere with the dipole trap- 
mechanism. However, the spin coupling, given by the 
term 



1m e 



S-B 



(44) 



in the Hamiltonian H' el , oscillates with the optical fre- 
quency, and then much faster than all the other rel- 
evant frequencies involved. The probability for spin 
transitions in first-order perturbation theory scales as 
(ge^s-Si) 2 / \fojOL) 2 ~ 10~ 15 (where the subscript L indi- 
cates that the magnetic field corresponds to the trapping 
laser). Second-order processes have associated rates of 
the order of ~ I0 -15 wl, and thus very slow compared to 
the other dynamical phenomena under study. 

The relevance of the remaining relativistic corrections 
to the coupling between the laser field and the ion can 
be estimated by plugging orders of magnitude and using 
parity and resonance arguments, as has been done in the 
previous subsections. The largest term comes from the 
presence of the electric field of the laser in the spin-orbit 
coupling in the third line of eq. (|30|) . which gives rise to 
an interaction of the form: 



(2m e c) 



x p 



(45) 



This extra coupling can be added to the dipolar Hamil- 
tonian ([T]), causing a spin-dependent modification of the 
coupling constant of order at most 10" 5 . As was the 
case with the corrections found in the previous subsec- 
tions, this is a charge-independent effect, namely, it is 
present also in the trapping of neutral atoms. 



FINAL REMARKS 



These two couplings are in principle of the same order, 
but the first of them is off-resonant from the transitions 
it couples, so that the dominant contribution will be 
given by the second one. In a two-level description of 
the dipolar trapping, the latter term corresponds to a di- 
agonal operator which induces different energy shifts in 
the states \g), |e), thus modifying the frequency of the 
atomic transition. For the magnetic field of the trapping 
laser, the order of magnitude of this intensity-dependent 
shift is of about I Hz, which means nine orders of mag- 
nitude smaller than the shift associated to the potential 
depth. This effect is thus smaller than the ones found 
previously, and it is also independent of the total charge. 



We have studied the trapping of an ionized atom with 
one valence electron by means of weak optical fields. The 
different terms that couple the ion and the field have been 
analyzed, and the effects due to the net charge of the par- 
ticle have been identified. From the estimations carried 
out in this work, one can conclude that the couplings rel- 
evant for the dipolar trapping of a charged particle are 
to very good approximation equivalent to those present 
for neutral particles. The most important effect due to 
the charge monopole is associated with the direct cou- 
pling between the vector potential of the external field 
and the center-of-mass momentum of the charged par- 
ticle. For the field of the trapping laser, this term gives 
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rise to a micromotion-like effect with energy of order 10 -8 
with respect to the potential depth of the dipolar trap, 
as found in subsection IV (i) The net charge can also 



be affected by other external electric fields, such as the 
one used in [l2| to confine the particles in the axial direc- 
tion. This implies the need to carefully suppress unde- 
sired additional fields, specially close to resonance with 
the center-of-mass motion; however, only dc-components 
of the stray fields get compensated up to now. Black- 
body radiation at room temperature, nevertheless, does 
not couple strongly enough to the charge so as to produce 
significant heating in the time scales of the problem. 

One interesting aspect that remains open is connected 
with the minimum temperatures achievable in dipole 
traps when ions and atoms are trapped together. In- 
deed, one of the motivations for the optical trapping of 
ions is the possibility to use the same setup to confine 
both neutral and charged particles forming composite 
systems. Such hybrid systems have already been exper- 
imentally realized with ions in Paul traps, allowing for 
the observation of cold ion-atom collisions [lj] and sym- 
pathetic cooling of an ion by a Bose-Einstein condensate 
15]. However, the micromotion of the ions in the Paul 
trap represents a limitation to the range of energies that 
can be explored [TU, [lTj . The possibility of optical trap- 
ping of ions is then an attractive alternative, since the 
micromotion of an ion in a dipolar trap is much smaller 
than in rf traps. 

Nevertheless, even in purely optical traps, only the 
charged particles will be subject to the extra driving due 
to the coupling of the monopole with the field. The effect 
of this coupling has been studied in Subsection IV (i) and 



the kinetic energies associated to it were estimated to be 
of about (QA) 2 /(2M). This phenomenon can thus give 
rise to atom-ion collisions in a similar way as in the case 
of Paul-trap setups, but with an energy scale which is 



smaller by several orders of magnitude: Micromotion in 
radiofrequency traps has associated energy scales of the 
order of ks times a few /iK (for very cold ions and good 
micromotion compensation), whereas the driven motion 
in the dipolar trap for the setup of [l2[ can be estimated 
to have a kinetic energy smaller than ks X 1 nK. We note 
that for red-detuned trapping beams both the kinetic 
energy associated to this micromotion and the trapping 
depth are proportional to the laser intensity; therefore, if 
the laser intensity is increased to achieve larger trapping 
depths the energy associated to the driven oscillation will 
increase by the same amount. However, for blue-detuned 
trapping the particles are confined around the nodes of 
the field, and thus this micromotion-like effect can be fur- 
ther suppressed. These considerations indicate that laser 
trapping, for instance in optical lattices, is a promising 
approach for realizing hybrid ultracold atom-ion systems 

The realization of composite charged-neutral systems 
can allow for the study of impurities in a bath [36[ and ul- 
tracold atom- ion collisions [13] ■ The possibility to create 
optical potentials for ions has also been proposed for the 
implementation of Frenkel-Kontorova models [38]. We 
expect our work to provide a useful reference for an ac- 
curate description of the corresponding dynamics. 
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